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摘 要 : 表 奢 履 盖 型 冰川 在 我 国 西部 分 布 广泛 ,由 于 该 类 型 冰川 消融 区 域 不 同 程度 的 被 岩石 碎 丑 所 履 盖 ,消融 状 
况 与 非 表 厂 覆盖 型 冰川 有 较 大 差异 ,因此 ,开展 表 厂 履 盖 型 冰川 的 消融 模拟 研究 至 关 重 要 。 本 文 以 冰 面 气象 数据 


为 驱动 ,使 用 表 磺 覆盖 型 冰川 能 量 平衡 模型 对 天 山 托 木 尔 ! 


传导 过 程 和 能 量 平衡 方程 ,模型 计算 了 表 磺 表面 温度 以 及 表 磺 层 内 部 温度 ,并 通过 表 磺 内 部 温度 估算 下 鹤 冰 的 消 
融 量 。 结 果 表 明 : 2008 年 夏季 期 间 模 型 模拟 的 消融 量 为 0.39 m w.e., 与 消融 花 杆 数据 进行 验证 取得 了 较 高 的 模拟 


条 青 冰 滩 72 号 冰川 表 磺 区 进行 能 量 和 消融 模拟 。 基 于 热 


精度 (R=0.92, RMSE=+0.03 m w.e.) , 表 磺 表面 温度 和 内 部 10 cm 深 处 温度 的 模拟 值 也 与 实测 数据 拟 合 较 好 ( 尼 分 别 


为 0.91 和 0.60)。 在 表 太 区 的 能 量 交换 过 程 中 , 净 短 波 辐射 是 


的 能 量 收入 项 , 感 热 通 量 是 最 大 的 能 量 支 出 项 


(49.796) ,其 次 分 别 为 传导 热 通 量 (消融 耗 热 )(25.8%), 净 长 波 辐射 (19.8% ) 和 潜 热 通 量 (4.6%) ,降水 热量 不 足 1%。 
云 量 对 表 研 区 的 气象 和 能 量 特征 有 着 显著 的 影响 , 阴 天 条 件 下 表 奔 区 的 人 射 短波 辐射 峰值 从 晴天 的 854 W «m? KE 


E 


28587 Wm”, 下 行 长 波 辐射 和 相对 湿度 增加 ,平均 消融 量 
析 表 明 ,模拟 的 消融 量 对 导热 系数 的 变化 最 为 敏感 ,反照 率 和 不 


量 比 晴天 减少 了 12%。 此 外 ,对 表 研 关键 参数 的 敏感 性 分 


1 表面 粗糙 度 的 变化 量 同样 不 可 忽视 。 


关键 词 : 青 冰 滩 72 号 冰川 ; 天 山 托 木 尔 峰 ; doi. 能 量 平衡 ; 消融 模拟 


20 世 纪 以 来 ,全 球 变 暖 导致 的 冰川 加 速 退 缩 已 


能 量 传导 至 冰 面 时 消耗 过 多 ,进而 抑制 下 获 冰 的 消 


引起 了 广泛 关注 ,山地 冰川 作为 中 亚 干 旱 内 陆 区 域 
的 重要 淡水 来 源 , 其 消融 变化 与 融 水 径流 显著 影响 
下 游 的 农业 和 生产 生活 用 水 。 中 国境 内 表 碘 覆 
六 型 冰川 共有 1723 条 ,总 面积 为 12974 km , 表 磺 所 
占 面积 为 1494 km*(11.5%); 尤 其 是 在 我 国 天 山 , 表 
磺 覆 盖 型 冰川 数量 最 多 ,面积 也 最 大 ,所 占 全 部 冰 
川 面积 比例 达 54% , 表 研 覆盖 区 域 约 481 km’, zi 


融 。 因 此 ,研究 表 矿 的 物理 特征 (如 厚度 .导热 系 
数 .表面 粗糙 度 ) 和 能 量 平衡 过 程 对 揭示 表 矿 履 盖 
型 冰川 的 消融 状况 尤为 重要 。 

表 厂 覆盖 型 冰川 能 量 平衡 模拟 一 直 被 学 者 们 
所 关注 。Nakawo 等 5 基于 野外 表 夺 观测 数据 计算 
了 表 磺 下 覆 冰 的 消融 量 ,得 到 了 与 实测 消融 量 较为 
一 致 的 模拟 结果 ,并 给 出 了 通过 表 矿 热 阻 估算 消融 


区 域 占 表 厂 覆盖 型 冰川 面积 的 12.2%"。 表 矿 的 存 
在 使 得 冰 / 雪 表面 与 大 气 的 能 量 交换 变 为 表 厂 表面 
与 大 气 的 能 量 交 换 ,其 消融 特征 显著 区 别 于 裸 冰 或 
积 雪 表 面 。 大 量 研究 表明 =" , 表 厂 厚度 小 于 某 一 阔 
值 时 ,由 于 表 厂 表面 的 反照 率 小 于 雪 冰 表面 致使 表 
磺 吸收 更 多 的 太阳 辐射 ,从 而 加 速 表 磺 下 覆 冰 的 消 
融 ; 当 表 柄 厚度 大 于 这 一 闵 值 时 , 表 磺 表面 鱼 余 的 
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量 的 方法 ;Nicholson 等 "使 用 基于 表 研 内 热传导 模 
块 的 能 量 平衡 模型 ,结合 日 尺度 气象 数据 作为 模型 
驱动 数据 ,应 用 在 阿尔 插 斯 山 和 斯 瓦尔 巴特 群岛 两 

条 冰川 ,结果 表明 模拟 的 下 禾 冰 消融 量 较 好 地 反映 
了 实测 消融 量 。 此 后 , 表 磺 覆盖 型 冰川 的 模拟 研究 
基于 野外 实测 数据 ,开展 了 表 矿 层 垂 直 结 构 内 的 水 
汽 和 热传导 等 更 为 细致 的 研究 。 其 中 ,Rounce 
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等 "通过 对 表 奢 导热 系数 .反照 率 和 表面 粗糙 度 进 
行 了 野外 测量 ,在 此 基础 之 上 模拟 了 表 磺 表面 温度 
和 下 和 覆 冰 消融 量 ,分 析 了 模型 参数 的 不 确定 性 对 模 
拟 结 果 的 影响 程度 。 国 内 有 关 表 磺 覆 盖 型 冰川 的 
模拟 研究 起 步 相 对 较 晚 ,Han Se SEP a BUG E 
理论 构建 了 适用 于 天 山 科 其 喀 尔 冰川 表 磺 状况 的 
能 量 平衡 模型 ;Yang 等 "利用 表 磺 能 量 平衡 模型 分 
析 了 藏 东南 号 隆 拉 冰川 的 气象 和 辐射 分 量 特征 ,与 
非 表 磺 和 覆盖 型 的 帕 隆 4 号 冰川 相 比 , 嘎 隆 拉 冰 川 吸 
收 了 更 多 的 太阳 辐射 和 下 行 长 波 辐射 ,从 而 增加 了 
上 行 长 波 辐射 和 潮流 热 通 量 的 亏损 进而 导致 消融 
量 增多 。 此 外 ,Zhang 等 ”构建 了 基于 表 矿 热 阻 系 
数 的 能 量 平衡 模型 ,进一步 发 展 了 表 磺 厚度 估算 模 
型 ,并 在 藏 东南 表 奢 覆盖 型 冰川 的 分 布 区 域 得 到 了 
应 用 。 基 于 上 述 研 究 发 现 , 表 厂 履 盖 型 冰川 的 模拟 
研究 主要 集中 在 阿尔 插 斯 山脉 和 青藏 高 原 地 区 ,天 
山 以 及 托 木 尔 峰 地 区 表 研 覆盖 型 冰川 的 消融 研究 
较 少 ,限制 了 我 们 对 该 区 域 冰川 消融 速率 和 冰川 融 
水 资源 变化 的 认识 。 因 此 ,本 文 从 冰川 能 量 平衡 过 
程 出 发 ,对 青 冰 滩 72 号 冰川 进行 消融 模拟 研究 , 验 
证 能 量 平衡 模型 和 表 磺 属性 参数 在 实际 模拟 中 的 
适应 性 ,完善 对 该 区 域 表 研 覆盖 型 冰川 消融 机 制 的 
认识 。 

天 山 托 木 尔 峰 地 区 作为 亚洲 内 陆 最 大 的 冰川 
发 育 区 ,该 区 蕴藏 着 丰富 的 冰雪 水 资源 , 托 木 尔 峰 
山区 平均 每 年 产 流 约 63.4xl0'm ,其 中 约 50% 为 冰 
crit z^ ,对 区 域 水 资源 和 人 民生 产生 活 至 关 重 
要 ,同时 对 该 区 域 冰川 的 现状 和 消融 研究 娘 终 是 国 
内 外 冰川 学 的 研究 热点 "。 鉴 于 此 ,本 文选 择 天 山 
托 木 尔 峰 青 冰 滩 72 号 冰川 作为 研究 对 象 ,基于 实测 
气象 数据 ,结合 野外 表 奢 温度 和 消融 量 观测 ,使 用 
表 矿 覆盖 型 冰川 能 量 平 衔 模 型 开展 表 厂 区 冰川 消 
融 模拟 研究 ,对 表 磺 区 的 能 量 平 衡 特 征 和 热传导 状 
况 进行 分 析 , 以 期 揭示 表 夺 覆 盖 对 冰川 能 量 收 文 和 
消融 过 程 的 影响 。 


1 研究 区 概况 


托 木 尔 峰 地 区 位 于 中 国境 内 天 山 山脉 的 最 西 
端 ,是 亚洲 内 陆 最 大 的 冰川 发 育 区 之 一 。 据 中 国 冰 
川 编目 资料 , 托 木 尔 峰 地 区 在 我 国境 内 共有 冰川 
1858 4& ,总 面积 达 4195 km? ,是 我 国 新 疆 阿 克 苏 地 


区 和 伊犁 地 区 重要 的 水 资源 ,也 是 塔里木 盆地 北 侧 
塔里木 河 主要 支流 的 源 区 。 托 木 尔 峰 地 区 为 大 陆 
性 干旱 气候 ,冰川 区 雪线 附近 年 平均 气温 为 -11~ 
-7% ,降水 主要 来 自 大 西洋 和 北冰洋 的 暖 湿 气流 ， 
RAA 5—9 H 

HYKIE 72 c vk JI] (41°45.51'N, 79*54.43' E) 发 
DEF Bey oe aN TAT I Sc Ut EES Hie TAY EAE DK Sa, Hu FE 
木 尔 峰 南 端 ,是 一 条 朝 南 的 山谷 - 冰 斗 冰川 ,由 降雪 
和 冰 / 雪 月 补给 (图 1)。 青 冰 滩 72 号 冰川 面积 》 
5.61 km , 冰 舌 末端 至 冰川 顶部 的 海拔 介 于 3560~ 
5986 m 之 间 , 冰 川 长 度 约 为 7.4 km, 冰 川 存在 冰 崖 和 
DK TET CU. Petia EEE Ait Y d D Ok TAS Shin Rl 
T3 A), Ze fi FE BE AT FAY 0.87 km? , 占 消 融 区 面积 的 
52%1'””1。 表 磺 厚 度 整 体 随 海拔 升 高 而 减少 ,厚度 
最 大 处 超过 1 m, 冰 舌 东 侧 和 西 侧 表 研 的 平均 厚度 
分 别 为 7cm 和 16 em, H PIMA Rm he A FRA EAR 
侧 。72 号 冰川 的 表 研 临界 厚度 约 为 4cm, 超 过 该 厚 
度 的 表 矿 区 消融 量 随 厚度 增加 而 减少 ,其 面积 约 为 
0.66 km’, 占 表 奢 总 面积 的 76% ; 表 厂 厚度 小 于 4 cm 
的 区 域 消融 量 随 表 矿 厚度 的 增加 而 变 大 ,面积 ， 
0.21 km22 与 天 山 山脉 其 他 冰川 相 比 ,该 冰川 运 
动 速度 较 快 ,年 均 水 平 运动 速度 为 47 ma 


2 数据 与 方法 


气象 数据 

本 文 的 气象 数据 来 源 于 在 青 冰 滩 72 号 冰川 表 
磺 区 架设 的 自动 气象 站 (AWS ) ,型 号 为 Davis Van- 
tage Pro2 Plus。 该 气象 站 位 于 冰川 冰 舌 东 侧 海拔 
3950 m lb ,该 处 的 表 矿 厚度 为 12 cm, 主要 由 灰色 和 
白色 花岗岩 及 碎 导 组 成 (图 1)。 观 测 的 气象 要 素 包 
括 气 温 、 相 对 湿度 、 风 速 、 降 水 、 入 射 /反射 短波 辐射 
等 ( 表 1) ,采集 时 间 间 隔 为 1h。 本 研究 使 用 的 气象 数 
据 时 段 为 2008 年 7 月 30 日 14:00 一 8 月 29 H 22:00, 
由 于 复杂 冰 面 气象 环境 的 影响 ,观测 期 间 自 动 气象 
站 的 数据 记录 在 2008 年 8 月 6 日 17:00 一 7 H 12:00 
和 8 月 8 日 14:00 一 16:00 期 间 存 在 一 定 的 缺失 。 需 
要 说 明 的 是 ,由 于 缺失 的 时 间 段 内 未 能 进行 消融 模 
拟 , 所 以 数据 缺失 时 间 段 的 消融 量 由 平均 消融 量 代 
TR. ERF, HH T fm a TE 2008 4 8 H 28 H 
12:00—8 H 29 日 22:00 之 间 无 记录 ,因此 ,模拟 的 
结束 时 刻 提 前 了 34 个 小 时 。 
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TE: (a) 研究 区 地 理 位 置 ;(b) 青 冰 滩 72 号 冰川 表面 地 形 ;(c) 自动 气象 站 ;(d) 消融 花 杆 。 
1 研究 区 概况 示意 图 
Fig. 1 Map of study area 
表 1 自动 气象 站 气象 传感器 和 指标 
Tab. 1 Meteorological sensors and their indicators of the automatic weather station 
测量 项 目 气温 /% 相对 湿度 /% 风速 /ms ) 降水 /mm 入 射 /反射 短波 辐射 /((W.m”) 
精度 +0.5 *C +3% <4% +4% +5% 
架设 高 度 2 m 2 m 2 m 0.7 m 1.5m 


2.2 ERAS 气 象 再 分 析 数 据 

ERA5 是 ECMWF 对 过 去 全 球 气候 和 天 气 的 第 
五 代 再 分 析 产 品 , 其 水 平分 辨 率 为 0.25°x0.25° ,时 
间 分 辩 率 为 1h。 本 文选 取 了 研究 期 内 ERA5 云 量 
因子 数据 分 析 不 同 天 气 状况 下 气象 和 能 量 条 件 的 
变化 。 不 同 的 天 气 状况 划分 参考 Van Den Broeke 
等 ”基于 云 量 因子 的 划分 方法 ,ERA5 的 云 量 因 
子 <0.3 时 划分 为 晴天 , > 0.7 时 划分 为 阴 天 ,不 考虑 
降水 的 影响 ,划分 后 通过 逐 小 时 的 能 量 分量 时 间 序 
列 变 化 进行 分 析 。 数 据 来 源 于 哥 白 尼 气 候 变 化 服 
务 (Copernicus Climate Change Service) 数据 平台 
(https://eds.climate.copernicus.eu) « 
2.3 slim E RR 

Be tot ia BE Ud fe e oi FL sc AL YK JII BE te F pp 
拟 的 关键 参数 和 模拟 结果 的 重要 验证 资料 。 研 究 
人 员 于 消融 期 间 在 冰川 海拔 3950 m b B] Roi DC 
了 一 个 2 症 的 温度 孔 。 该 温度 孔 处 表 犬 层 厚 度 为 


13 cm ,在 表 厂 层 内 共 布 设 了 2 个 温度 探头 。 其 中 一 
个 探头 置 于 表 奔 表面 以 下 约 1cm 人 处 ,用 于 测量 表 碘 
表面 的 温度 。 由 于 表面 温度 的 观测 容易 受到 太阳 
辐射 的 影响 ,因此 使 用 表 磺 表面 以 下 1 cm 处 的 温度 
来 近似 代替 表面 温度 ,这 与 Rounce 等 "测量 喜 马 拉 
雅 山 脉 tmja-Lhotse Shar 冰川 表 奔 表面 处 温度 的 方 
法 相 一 致 。 另 一 个 温度 探头 布设 在 表 磺 表面 以 下 
10 em 处 ,用 于 测量 表 厂 层 内 部 的 温度 。 表 磺 以 下 
冰 体 内 ,每 隔 20 em 布设 一 个 温度 探头 测量 冰 温 。 
测 温 探头 均 为 热 敏 电阻 ,由 中 国 科 学 院 西北 生态 环 
境 资源 研究 院 ( 原 中 国 科学 院 寒 区 旱 区 环境 与 工程 
人 研究 所 ) 研 制 , 精 度 为 :0.05 CCC。 野外 试验 中 ,人 工 使 
用 万 用 表 测 量 探头 电阻 率 ,进而 计算 表 磺 层 和 冰 体 
内 部 的 温度 。 共 开展 了 6 次 观测 ,分 别 为 7 月 31 日 、 
8 月 1 日 8 月 3 日 8 月 4 日 8 月 7 日 和 8 月 8 日 。 
2.4 冰川 实测 消融 量 数据 

研究 人 员 在 自动 气象 站 附近 表 厂 区 布设 了 1 根 
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消融 花 杆 (图 lc) , 表 厂 厚度 为 12 em。 在 2008 年 7 
月 29 日 一 9 月 1 日 期 间 , 开 展 了 26 次 详细 观测 。 具 
体 的 观测 内 容 为 花 杆 至 冰川 表面 的 垂直 高 度 。 消 
融 量 的 计算 公式 为 : 


A 2p, ^m (1) 
式 中 ;4 为 消融 量 (m w.e.);pi; 为 冰川 冰 的 平均 密度 
(900 kg- m™°); Am 为 观测 期 间 花 杆 至 冰 面 的 垂直 高 
度 变 化 (m)。 最 终 得 到 研究 期 内 花 杆 的 消融 量 值 ,本 
研究 使 用 该 数据 作为 冰川 消融 量 模 拟 的 验证 资料 。 
2.5 表 磺 区 能 量 平 衡 方 程 
本 文 使 用 表 磺 覆盖 型 冰川 能 量 平衡 模型 (De- 
bris Energy-Balance model ,简称 DEB 模 型 ) 进 行 单 点 
尺度 冰川 能 量 和 消融 模拟 ,该 模型 开发 较为 成 熟 ， 
物理 意义 明确 ,完整 描述 了 表 磺 区 冰川 的 能 量 交换 
过 程 ,与 度 日 模型 相 比 能 获取 能 量 分量 对 消融 过 程 
的 影响 , 且 作 为 一 款 开 源 模型 在 阿尔 插 斯 山脉 
Miage 冰川 和 藏 东南 嘎 隆 拉 冰 川 取得 了 很 好 的 模拟 
AR, HP , 表 夸 覆盖 型 冰川 能 量 平 衔 模 型 主 
要 基于 如 下 能 量 平衡 等 式 : 
S+L+H+LE+P+G+M=0 (2) 
式 中 :5 为 净 短 波 辐射 ;L 为 净 长 波 辐 射 ; 甩 为 感 热 通 
量 ;LE 为 潜 热 通 量 ;P 为 降水 带 来 的 热 通 量 ;G 为 传 
导热 通 量 , 即 能 量 从 表 硬 表面 传输 至 冰 面 所 需 的 损 
耗 ;WM 为 冰川 消融 能 量 ;单位 均 为 Wm ”*。 在 本 文 的 
模型 中 , 表 奔 表面 温度 和 上 行 长 波 辐射 整体 带 入 能 
量 平衡 方程 求解 , 表 研 表 面 温度 7 不 仅 是 计算 各 辐 
出 通 量 的 输入 项 ,同时 其 数值 也 取决 于 辐射 通 量 的 
大 小 。 除 太阳 短波 辐射 和 下 行 长 波 辐 射 外 ,其 他 能 
量 通 量 均 可 以 表示 为 表 磺 表面 温度 7 的 函数 。 该 模 
型 使 用 Crank-Nicolson 数值 迭代 方案 ,计算 时 间 步 长 
为 1h,7 既 是 能 量 平衡 模型 中 的 未 知 数 , 又 是 计算 能 
量 项 的 所 需 参 数 , 通 过 将 表 磺 分 为 数 层 EC 
层 表 磺 的 温度 ,从 而 得 出 冰 面 的 热 通 量 "””。 由 于 表 
磺 之 间 的 缝 除 中 存在 大 量 的 空气 和 水 汽 ,太阳 辐射 
和 淇 流 热 通 量 在 表 厂 内 发 生 复杂 的 能 量 交换 ,能量 
平衡 模型 根据 表 磺 内 的 温度 变化 表示 能 量 传递 过 
程 。 受 太阳 辐射 和 长 波 辆 射 的 影响 , 表 犬 表面 温度 
时 刻 发 生变 化 ,昼夜 温差 大 。 但 根据 表 夸 温度 观测 
和 自动 气象 站 气温 数据 , 表 磺 层 温 度 在 消融 期 内 基 
本 都 为 正 温 ,因而 可 假定 表态 底部 与 冰 体 交界 处 的 
温度 一 直 处 于 融 点 (0%)。 与 此 同时 ,通过 对 表 矿 
不 同 厚 度 处 的 温度 进行 观测 , 表 厂 内 的 温度 曲线 呈 


现 出 非 线性 变化 , 异 夜 的 能 量 传递 方向 也 可 能 发 生 
改变 。 因 此 ,使 用 非 定 态 热传导 方程 描述 表态 内 的 
温度 变化 : 


QST _ oT (3) 
az ðt 
a= (4) 
Pala 
式 中 :ao 为 热 扩 散 系数 ,表示 非 定 态 热传导 过 程 中 物 


体内 部 温度 趋 于 均匀 的 能 力 ;7 为 某 一 表 厂 层 的 温 
度 (K);z 为 每 一 表 奢 层 的 深度 (由 ;为 时 间 (h 因 为 表 
fit FINA BW mm! K’), pa Fc BN Be tk AY BE 
(1496 kg'm”) 和 热 容 量 (948 J kg K’) o 
2.5. 短波 辐射 ”地 表 的 短波 辐射 来 源 于 太阳 辐 
射 , 同 时 也 是 冰川 表面 的 主要 能 量 来 源 。 云 量 对 入 
射 短波 辐射 有 重要 的 影响 , 表 磺 反照 率 很 大 程度 上 
也 决定 了 出 射 短波 辐射 的 大 小 。 本 文 使 用 自动 气 
象 站 的 入 射 /出 射 短 波 辐射 数据 作为 模型 输入 。 
2.5.2 长 波 辐射 ”长 波 辐射 包括 地 表 发 射 的 上 行 长 
波 辐射 和 大 气 发 射 的 下 行 长 波 辐射 ,由 于 自动 气象 
站 处 没有 长 波 辐射 的 实测 数据 ,因此 使 用 斯 特 藩 - 
玻 尔 兹 曼 定律 计算 下 行 和 上 行 长 波 辐 射 ””*。 下 行 
长 波 辐射 应 的 计算 公式 为 : 

Losegsl. (5) 
式 中 :o 为 玻 尔 兹 曼 常 数 [5.67x10*W: m K]; T, 
为 大 气温 度 (K);es 为 有 效 辐射 率 , 即 考虑 云层 覆盖 条 
件 下 的 大 气 辐射 率 ,采用 Idso” 和 Unsworth 等 ”的 公 
式 进行 计算 : 


1500 


e 20.74 5.95x 107e*e " (6) 
£ a =8 (1 - 0.84n) + 0.84n (7) 
SU Lei NEZ ED ses D TARE TE TK IR (Pa) s n W 
云 量 , 来 自 于 ERAS 再 分 析 数 据 。 

上 行 长 波 辐射 L, 的 计算 同样 采用 黑体 辐射 定 

律 ,其 不 确定 性 主要 取决 于 模拟 的 表 奔 表面 温度 : 
L,76,0T; (8) 

式 中 :se 为 表 矿 的 比 辐射 率 ; 了 为 表 厂 表 面 温度 (四 )。 
2.5.3 Wü jLihi8 X 低层 大 气 与 下 热 面 之 间 因 温度 
或 湿度 梯度 产生 的 能 量 交 换 称 为 满 流 热 ,包括 感 热 
通 量 和 潜 热 通 量 。 注 流 热 通 量 是 冰川 能 量 平衡 的 
重要 组 成 部 分 。 通 过 计算 整体 理 查 森 数 和 判断 冰 
面 稳定 状态 可 以 较 好 的 模拟 洁 流 热 通 量 , 因 此 ,本 
文 使 用 整体 空气 动力 学 方法 计算 感 热 通 量 ,计算 所 
需 的 表 厂 表面 温度 数据 同样 是 模拟 值 ,其 基本 思路 
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SAAB ia Pt 2388 ec 5 t BE RE ,湿度 梯度 和 风速 梯 
度 存在 相关 关系 。 感 热 通 量 是 近 地 面 气温 与 下 
热 面 温度 之 间 的 差异 引起 的 热量 交换 。 潜 热 通 量 
则 是 近 地 面 大 气 与 下 热 面 之 间 水 汽 相 变 产生 的 能 
量 交 换 , 包 括 燕 发 和 凝结 。 感 热 通 量 及 和 潜 热 通 量 
LE 的 计算 公式 为 : 


T =T = 
H-p,c ku x : (5,,) : (9) 
(I z Jz] 
20m Zor 
LE=p,L, ku ai p A (10) 


(m) In 
Zo m 204 


StH sg Ae NE (m * s’); T, 29 CAR X 8 C 
7 与 表 磺 表面 气温 7 的 平均 值 。 
2.5.4 降水 热量 ”降水 与 表 厂 之 间 的 热量 交换 已 采 
用 Reid 等 ”提供 的 方法 计算 : 

P =p,c,r(T, 一 T.) (14) 
SUF ip, HK BY PE BE (kg * m7); e, WIK B EG 3 (4.18 
10° 本 kg C 7); 为 降雨 速率 (ms1') ,了 为 降水 温度 
() ,由 于 缺乏 实地 观测 ,使 用 气温 已 代 替 降 水 温度 。 
2.5.5 传导 热 通 量 ”一般 将 物体 或 系统 内 的 温度 差 
产生 的 热量 传递 现象 称 为 热传导 。 导 热 系数 上 是 计 
算 传 导热 通 量 的 关键 参数 , 它 用 来 描述 物体 传导 热 
量 的 固有 能 力 ,热量 传递 的 速率 取决 于 温度 梯度 的 


式 中 :p, 为 当地 空气 密度 (kg' mic 为 空气 热 容量 
("kg + K); LIRR * kg '), 近 地 面 气 温 为 
283 KK 时 其 数值 为 2.476x10J'kg ';hki 为 卡门 常数 
(Von Karman constant) ,是 假设 混合 长 和 速度 廓 线 的 
关系 所 引入 的 经 验 系 数 ,实验 室 和 近 地 面 大 气 测定 
值 在 0.35~0.43 之 间 ,本 文 取 值 为 0.4 记 为 近 地 面 风 
3E (m * s7); T, 2g Xr Hb. rfi ^C il (K); T. 29 Pe fiot de Th il SE 
(K); e Fl e, 4) 31 Ay VF Hh. i KAMK o e ELE KC TAL 
(Pa);z 为 自动 气象 站 气温 传 感 絮 距离 地 面 的 高 度 
(2 m) zo ,为 表 磺 表面 (动量 ) 的 粗糙 度 (m) ,定义 为 水 
平 风 速 为 0 处 距 地 物 表面 的 高 度 , 取 值 为 0.016 m^"; 
zo 和 zo 分 别 为 表 磺 表面 的 热量 和 水 汽 通 量 的 粗 烽 
度 , 二 者 取 值 与 a 相同 。@B,、8B, 和 ,分 别 为 动量 、 
热量 和 水 汽 的 无 量 纲 稳定 性 函数 。 在 计算 潜 热 通 
量 的 过 程 中 , 需 定 义 表 太 表面 的 相对 湿度 ,本 文 将 
自动 气象 站 处 测 得 的 降水 量 > 0.2 mm 时 的 表面 相 
对 湿度 定义 为 100% , 15 WU] Jy 096^; 
Brutsaert 等 使 用 理 查 森 数 尺 来 描述 动量 .热量 
和 水 汽 通 量 的 稳定 性 状态 ,本文 假 设 自动 气象 
站 处 风速 为 0 时 ,三 者 处 于 稳定 状态 ,否则 将 处 于 非 
稳定 状态 。 当 三 者 在 近 地 面 处 于 稳定 状态 时 , 理 查 
森 数 为 正 ,可 表示 为 : 
(BB) -(5,5) -(1-5R) (11) 
当 近 地 面 处 于 非 稳定 状态 时 , 理 查 森 数 为 负 ， 
此 时 可 表示 为 : 
(5,5,)'-(6,5) -(1-16R) ^ (12) 
P : HA ERR; fH] Brutsaert 描述 的 方法 计 
AE ,公式 如 下 : 


0.75 


gelT, 一 T )(z 一 fis) 


- (13) 
Tw 


R,= 


大 小 和 表 矿 的 热 特性 , 表 矿 的 导热 系数 与 其 颗粒 大 
小 .密度 . 表 矿 内 孔 院 湿度 和 温度 密切 相关 ,其 数值 
大 小 受到 周围 环境 的 影响 , 表 厂 内 的 温度 和 湿度 越 
大 ,导热 系数 越 大 ,传递 到 冰 面 的 能 量 越 多 。 模 型 
通过 对 表 矿 进行 等 距离 分 层 ,每 层 表 矿 的 厚度 为 
0.01 m, 使 用 传 里 叶 定 律 计算 表 矿 表面 和 表 矿 / 冰 面 
交界 处 的 传导 热 通 量 ,使 得 每 层 表 矿 的 温度 变化 可 
以 近似 表示 为 线性 ,具体 计算 公式 如 下 : 

C= (Sr) T Es -T, s) 
式 中 :Gi 为 表 奔 底部 传导 至 冰 面 的 热 通 量 (W:m”); 
TA(N-1) 和 分 别 为 第 NN 层 表 碘 温度 和 冰 体 熔点 温度 
( 即 冰 面 温度 ) 次 为 每 层 表 矿 的 实际 厚度 (m)。 
2.5.6 冰川 消融 ”能量 平 衡 模型 模拟 的 表 矿 以 下 冰 
融化 量 近似 为 冰川 消融 量 。 冰 川 的 消融 量 c 由 可 供 
vr iege MOR ARO: 


sM (16) 
7 pi, 


SUP p NIKE E(kg m?) LAMME * kg )。 
3 结果 与 分 析 


3.1 模型 验证 

表 太 表面 温度 反映 了 表 矿 与 近 地 面 大 气 之 间 
的 能 量 交 换 程度 , 表 磺 表面 的 能 量 主要 以 热传导 的 
形式 向 下 传递 ,基于 表 奢 表面 温度 和 表 矿 内 热传导 
过 程 可 估算 出 下 履 冰 的 消融 量 ,野外 观测 的 表 奢 温 
度数 据 是 模型 的 直接 验证 参数 。 本 文选 取 了 能 量 
平衡 模型 模拟 的 表面 温度 和 10 cm 深度 处 温度 同 实 
测 值 进行 验证 (由 于 8 月 6 一 8 日 有 两 段 时 间 气 象 数 
据 存在 缺失 ,数据 缺失 的 时 间 段 内 未 显示 能 量 平衡 
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的 模拟 值 ) ,结果 表明 , RR TAT TARE CA 2a) AL 
值 和 实测 值 具有 很 好 的 一 致 性 ( 尼 =0.91, RMSE= 
1.78 C) ,同时 表 奢 下 10 em 温度 (图 2b) 的 模拟 值 也 
与 实测 温度 值 拟 合 较 好 (R=0.60, RMSE=0.48 °C) 。 
尽管 表 厂 下 10 cm 深度 处 的 温度 模拟 精度 有 所 降 
低 , 但 从 表 厂 履 盖 型 冰川 的 能 量 交换 过 程 中 可 知 ， 
表 磺 表面 温度 通过 近 地 面 气象 条 件 计算 得 到 ,进而 
基于 表面 温度 计算 表 矿 内 不 同 深度 处 的 温度 , 
IE, Ri F 10 cm 深 的 温度 模拟 结果 包含 了 表面 温 
度 的 计算 误差 。 从 表 矿 温度 的 时 间 序 列 变化 中 可 
以 看 出 ,表面 温度 的 日 波动 较 大 ,白天 温度 一 般 超 
过 15 ,夜晚 降温 迅速 ,温度 较 低 ; 表 磺 下 10 em 深 
度 靠近 冰 面 ,日 温差 小 ,白天 温度 不 超过 5 % ,夜晚 
温度 接近 0 % 。 基 于 野外 观测 的 温度 数据 中 ,有 44d 
的 正午 表面 温度 超过 15 % ,10 em 深 处 的 实测 值 也 
相应 较 高 ;剩余 2 d 温 度 低 于 10 % ,10 cm 深 处 的 温 
度 同样 较 低 。 通 过 对 气温 和 表 奔 表面 温度 进行 相 
关 性 检验 ,发 现 二 者 存在 较 好 的 相关 性 (r=0.55,P< 
0.01) ,因此 , 表 磺 温度 变化 在 一 定 程 度 上 可 以 反映 
近 地 面 气温 的 时 间 变 化 。 

利用 夏季 模拟 的 消融 与 实测 的 消融 进行 对 比 ， 
如 图 3 所 示 , 结 果 显 示 ,7 月 30 日 一 8 月 28 日 期 间 ， 
模拟 的 冰川 累积 消融 量 为 0.39 m w.e., 实测 的 累积 
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图 2 模拟 与 实测 的 表 厂 表面 温度 和 表 磺 下 10 cm 处 温度 
Fig.2 Simulated and measured debris temperature at Surface 
and depth of 10 cm 


消融 量 为 0.40 m we.。 虽 然 研 究 期 内 ,不 同时 间 段 
的 累积 消融 量 模拟 值 普遍 较 实 测 值 要 小 ,但 整体 来 
看 ,模拟 值 和 实测 值 两 者 之 间 具 有 较 好 的 一 致 性 ， 
决定 系数 ( 尼 ) 为 0.92, 均 方 根 误差 (RMSE ) X +0.03 
m w.e.。 两 者 之 间 存 在 差异 的 原因 可 能 主要 有 三 个 
方面 : 其 一 ,本 文中 表 厂 属性 参数 并 未 开展 实地 观 
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累积 消融 量 /(m w.e.) 


de 
X 
d 


时 间 / 月 -日 


SR 


图 3 模拟 与 实测 的 累积 消融 量 


Fig.3 Simulated and measured glacier melt 


测 ,是 通过 率 定 得 到 的 ;其 二 ,模型 没有 涉及 表 矿 内 
水 汽 的 蒸发 .再 冻结 过 程 , 模 拟 结果 与 真实 冰 面 能 
量 收 支 状况 存在 一 定 差距 ;其 三 ,该 条 冰川 运动 速 
度 较 快 ,存在 较 强 的 运动 补给 。 模 拟 与 实测 的 累积 
消融 量 在 8 月 15 一 19 日 间 均 呈现 出 较为 平缓 的 趋 
势 , 通 过 对 比 气温 数据 发 现 这 一 时 段 气温 整体 偏 
低 , 最 低 气 温 为 -3.9 C, RH ,此 时 气象 条 件 不 利于 
消融 ,以 致 消融 亏损 减少 。 尽 管 消融 的 观测 时 间 只 
有 一 个 月 左右 ,但 取得 的 观测 结果 精度 较 高 。 消 融 
观测 采用 的 是 花 杆 法 ,研究 人 员 每 隔 1~2 d 观 测 , 充 
分 记录 了 短 时 间 尺 度 内 消融 的 详细 变化 情况 ,为 开 
展 冰川 能 量 消融 模拟 提供 了 良好 的 实测 验证 数据 。 
3.2 RAK AS ESE 

Fa UK DE 72 7S vK)1] E AE dé DX BE tt REGIE WI 4 
和 图 5 所 示 。 在 亲 面 的 能 量 交换 过 程 中 , 
射 是 唯一 的 能 量 收入 项 。 感 热 通 量 是 最 大 的 能 量 
支出 项 (49.7%) ,其 次 分 别 为 传导 热 通 量 ( 消 融 耗 
热 ) (25.8% ) , f K UE fa BY (19.896 ) 和 潜 热 通 量 
(4.6% ) ,降水 热量 不 足 1% ,可 忽略 其 对 冰 面 能 量 交 
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图 4 表 研 表面 辐射 分 量 的 日 平均 变化 


Fig. 4 Mean diurnal variation of debris surface heat fluxes 


components 
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1250 上 一 一 净 短 波 辐射 ”一 一 传导 热 通 量 一 潜 热 通 量 
一 一 净 长 波 辐射 ”一 一 感 热 通 量 


能 量 值 (W'm 3) 
N 
3 


时 间 / 月 -日 


图 5 表 磺 表面 能 量 分 量 的 逐 小 时 变化 


Fig. 5 Hourly variation of debris surface heat fluxes 


components 


换 的 影响 。 

观测 期 内 净 短 波 辐射 表现 出 明显 的 日 变化 , 单 
日 峰值 可 超过 1000 Wm ?, 也 可 小 于 400 W:m’, iX 
可 能 受到 云 量 的 影响 ,具体 表现 为 云层 的 遮挡 大 幅 
减 小 了 入 射 短波 辐射 透 过 大 气 到 达 冰 川 表面 ,使 得 
净 短 波 辐射 曲线 出 现 扰动 或 峰值 削减 。 有 昼夜 变化 
方面 ,从 日 出 开始 短波 辐射 持续 增 大 ,午后 可 超过 
600 Wm”, 随 之 短波 辐射 持续 减 小 直至 日 落 。 感 热 
通 量 是 冰 面 的 最 大 能 量 支出 项 ,这 是 由 于 太阳 辐射 
能 量 加 热 表 厂 使 表 厂 温度 显著 高 于 近 地 面 气温 而 
产生 热量 对 流 交 换 。 午 后 感 热 通 量 的 亏损 高 达 
-400 负 :m? ,值得 注意 的 是 ,由 于 夜间 表 厂 温度 较 
低 , 使 得 感 热 通 量 处 于 轻微 的 鱼 余 状态 。 传 导热 通 
` 净 长 波 辐射 和 潜 热 通 量 均 呈 现 出 热量 支出 状 
,其 数值 在 白昼 期 间 波 动 不 大 ,曲线 变化 较为 平 
稳 , 且 亏损 量 不 超过 -200 Wm?*。 观 测 期 间 隆 水 带 
来 的 热量 几乎 为 0, 所 以 ,降水 热量 对 冰 面 能 量 交 换 


c mm 


的 消融 变化 与 气温 变化 具有 高 度 的 一 致 性 (r=0.73， 
P «0.01) ,最 高 的 日 均 温 (7.8 % ) 同 样 出 现在 8 月 3 
日 ,日 均 温 最 低 (-0.1 SC ) 的 一 天 为 8 月 19 日 ,这 说 明 
了 表 矿 覆盖 型 冰川 的 消融 变化 特征 同样 受到 环境 
气象 条 件 尤其 是 气温 的 控制 。 
3.4 云 量 对 表 研 区 气象 和 能 量 条 件 的 影响 

为 进一步 研究 不 同 天 气 条 件 下 表 研 区 的 气象 
和 能 量 特征 ,本 文 分 别 探讨 了 晴天 和 阴 天 天 气 中 气 
象 要 素 的 日 均 变化 , 云 量 对 各 种 气象 要 素 都 有 不 同 
程度 的 影响 (图 6) , 云 量 对 和 人 射 短波 辐射 的 影响 主 
要 体现 在 削弱 其 峰值 ,晴天 正午 左右 人 射 短波 辐射 
的 数值 高 达 850 W * m^, BAA AN E 600 Wm?, 云 
层 对 下 行 长 波 辐射 的 作用 则 完全 相反 , 阴 天 的 下 行 
长 波 辐射 始终 大 于 晴天 ,但 二 者 在 一 天 之 中 的 变化 
趋势 保持 一 致 。 晴 天 气温 略 高 于 阴 天 ,变化 幅度 在 
1.5~6 % 之 间 。 相 对 湿度 并 无 明显 的 变化 规律 , 阴 天 
的 相对 湿度 水 平 处 于 较 高 水 平 ,平均 相对 湿度 为 
70% 左 右 ,晴天 上 午时 间 段 内 相对 湿度 较 小 , 仅 为 
50% ,随后 相对 湿度 逐渐 增 大 。 一 天 之 中 风速 的 变 
化 有 较为 明显 的 2 个 峰值 ,正午 左右 风速 最 小 , 快 到 
日 落 时 风速 最 大 ,明天 的 风速 变化 与 晴天 相 比 较为 
平缓 。 降 水 主要 发 生 在 下 午 , 阴 天 的 降水 量 明显 较 
大 ,最 高 可 达 0.7 mm.pm:。 由 此 可 以 看 出 , 冰 面 的 气 
象 条 件 在 云 量 的 影响 下 时 间 分 布 特征 和 数值 大 小 
都 发 生 了 相应 改变 ,也 会 对 表 矿 区 的 热量 和 能 量 交 
换 造 成 一 定 的 影响 。 

通过 云 量 区 分 表 磺 区 的 天 气 状 况 后 ,可 以 看 
出 ,不同 天 气 条 件 对 能 量 平衡 方程 的 各 辐射 通 量 也 
产生 了 相应 的 影响 (图 7 和 表 2)。 云 量 对 净 短 波 辐射 


过 程 中 的 影响 可 忽略 不 计 。 表 磺 表面 各 辐射 通 量 
的 日 均 变 化 特征 较为 一 致 , 即 在 夜间 能 量 波动 小 且 


通 量 的 影响 较 大 ,晴天 正午 时 刻 的 净 短 波 辐射 超过 
750 Wm”, 阴 天 的 净 短 波 辐 射 则 削减 为 530 W * m^? 


趋 近 于 0, 日 出 前 后 收入 和 支出 量 持续 增加 ,在 午后 
绝对 值 达到 最 大 ,随后 逐渐 减 小 。 
3.3 表 磺 区 消融 特征 

基于 青 冰 滩 72 号 冰川 表 磺 区 的 气象 数据 ,结合 
能 量 平衡 方程 计算 出 模拟 期 内 的 日 平均 消融 量 为 
0.55 mm w.e., 人 研究 期 内 冰川 始终 处 于 消融 状态 ,日 
均 最 大 消融 发 生 在 8 月 3 日 ,消融 量 为 0.88 mm w.e., 
日 均 最 小 消融 发 生 在 8 月 17 日 ,消融 量 为 0.23 mm 
w.e., 从 7 月 30 日 一 8 月 28 日 ,累积 消融 量 呈 现 出 持 
续 减缓 的 趋势 ,为 进一步 研究 冰川 消融 的 变化 特 
征 ,对 模拟 期 间 的 日 平均 气温 进行 计算 ,发现 每 日 


左右 ,与 研究 期 内 所 有 天 气 条 件 相 比 ,晴天 的 平均 
净 短 波 辐射 增加 了 11.9% , 阴 天 则 减少 了 20.5%。 同 
样 的 ,晴天 条 件 下 感 热 通 量 的 最 大 亏损 量 超 过 -500 
W:m”, 阴 天 则 削减 为 -330 Wm 左右 ,晴天 和 阴 天 
的 平均 感 热 通 量 分 别 为 109.04 W * m?,72.30 W:m 
云 量 的 增加 虽然 减弱 了 净 短 波 辐射 的 收入 和 感 热 
通 量 的 支出 ,但 二 者 的 时 间 序 列 变化 特征 却 几 乎 没 
有 发 生变 化 。 阴 天 条 件 下 净 长 波 辐射 潜 热 通 量 和 
传导 热 通 量 与 晴天 时 相 比 ,其 数值 均 减 小 了 50% 左 
右 , 睛 天 时 洪 热 通 量 的 平均 亏损 量 较 小 , 仅 为 -2.63 
Wm , 净 长 波 辐射 和 传导 热 通 量 的 平均 亏损 量 则 
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图 6 不 同 天 气 条 件 下 气象 要 素 的 变化 


Fig.6 Effects of different weather conditions on meteorological elements 
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图 7 不 同 天 气 条 件 对 辆 射 通 量 的 影响 


Fig.7 Effects of different weather conditions on radiation 


flux 


在 阴 天 条 件 下 明显 减少 。 由 于 观测 期 间 降 水 稀少 ， 
降水 热量 可 忽略 不 计 , 在 不 同 天 气 条 件 下 也 无 明显 
变化 。 晴 天 和 阴 天 的 平均 消融 量 分 别 为 0.58 mm w.e. 
和 0.51 mm w.e., 阴 天 的 消融 量 比 晴 天 减少 了 12%。 
总 体 而 言 , 云 量 的 多 少 对 能 量 通 量 的 数值 大 小 有 显 
著 影 响 ,但 未 改变 能 量 通 量 随时 间 的 变化 特征 , 阴 


天 时 净 短 波 辐射 的 收入 和 其 他 能 量 项 的 亏损 都 有 
所 减少 ,晴天 时 除 潜 热 通 量 亏损 减少 外 , 
射 的 收入 和 其 他 能 量 项 的 亏损 都 在 增加 。 


4 讨论 


4.1 敏感 性 分 析 

表 奢 属性 作为 表 厂 覆盖 型 冰川 能 量 -消融 模 拟 
的 重要 不 确定 性 因素 ,本 文通 过 对 模型 中 各 表 厂 属 
性 参数 进行 敏感 性 分 析 ,可 以 确定 影响 表 厂 有 覆盖 型 
冰川 能 量 -消融 模拟 的 关键 参数 ,为 后 续 进 一 步 研 
究 提供 可 靠 的 基础 。 在 敏感 性 测试 中 , 除 表 厂 反照 
率 (Qa=0.086) 通 过 气象 站 的 反射 /入 射 短波 辐射 计算 
得 出 ,其 他 表 厂 参数 的 取 值 均 通过 文献 128] 获 取 。 
本 研究 分 别 选 择 表 矿 导热 系数 () . 表 厂 表面 粗糙 度 
(zo) \ 表 研 比 辐射 率 (8) 和 表 奔 反照 率 (a)4 种 表 碘 属 
性 参数 用 于 评估 表 厂 属性 参数 对 消融 模拟 的 影响 
( 表 3)。 在 对 青 冰 滩 72 号 冰川 的 实地 观测 和 模拟 
中 , 表 磺 导热 系数 设 定 为 0.94 Wem K, KFAR 
数 增加 /减少 0.2 Wm'…K' 时 ,模拟 消融 量 平均 增 
加 /减少 17.8%; 当 导 热 系数 增加 /减少 0.4 Wm K 
时 ,消融 量 平 均 增 加 /减少 35.7%。 模 型 中 表 奔 表面 
粗糙 度 设 定 为 0.016m, 当 粗 糙 度 增加 /减少 0.01 m 
时 ,模拟 的 消融 量 减 少 4.6% 或 增加 9.3%。 表 矿 比 辐 


1039] 何 捷 等 : 托 木 尔 峰 青 冰 滩 72 号 冰川 表 磺 区 夏季 消融 模拟 研究 1603 
表 2 不 同 天 气 条 件 下 的 平均 能 量 平衡 参数 分 量 及 其 变化 
Tab. 2 Mean fuxes and variation under different weather conditions 
净 短 波 辐 射 净 长 波 辐射 感 热 通 量 潜 热 通 量 传导 热 通 量 降水 热量 消融 量 
W:m”) KW *m?) KW *m?) KW *m?) /(W*m?) /(W*m?) /(mm w.e.) 
所 有 天 气 185.23 —36.58 -92.08 —8.62 —-47.88 -0.06 0.56 
HA 207.34 —46.05 —109.04 -2.63 —49.61 0 0.58 
变化 幅度 11.996 -25.9% -18.4% 69.5% -3.6% - 3.6% 
DX 147.30 —25.06 —12.30 —6.58 —-43.16 -0.20 0.51 
变化 幅度 -20.596 31.596 21.596 23.7% 9.9% -70.0% -8.9% 


RI 表 厂 参数 的 敏感 性 分 析 


Tab. 3 Sensitivity analysis of debris parameters 


参数 原始 值 。 。 参数 变化 量 rg ice 
k=0.94 W*m K” 10.2 0.47 17.296 
40.4 0.53 33.196 
-0.2 0.33 —18.496 
-0.4 0.25 —38.396 
zo=0.016 m +0.01 0.38 —4.696 
+0.02 0.37 -7.7% 
-0.01 0.44 9.3% 
£20.94 +0.03 0.38 —4.206 
+0.06 0.37 —8.296 
—0.03 0.42 4.496 
—0.06 0.44 9.196 
a=0.086 +0.1 0.38 -5.6% 
+0.2 0.35 -11.5% 
-0.05 0.41 2.796 


射 率 的 取 值 为 0.94, f HRIH Ie 0.03 时 , 模 
拟 的 消融 量 平均 减少 /增加 4.3%; 当 比 辐射 紊 增加 / 
减少 0.06 时 ,消融 量 平均 减少 /增加 8.7% 。 表 矿 对 
可 见 光 的 吸收 能 量 较 强 ,反照 率 接近 于 0, 气 象 站 处 
测量 的 平均 表 厂 反照 率 为 0.086 ,其 数值 增加 0.1 或 
0.2 ,模拟 的 消融 量 将 减少 和.6% 或 11.5% 。 除 导热 系 
数 外 , 表 矿 的 粗糙 度 . 比 辐射 率 和 反照 率 与 冰川 的 
消融 均 呈 现 出 负 相 关 关 系 , 表 磺 导 热 系数 在 能 量 平 
衡 模 型 中 表现 为 最 敏感 的 参数 , 表 厂 反照 率 和 表面 
粗糙 度 具 有 一 定 的 敏感 性 ,但 二 者 在 表 磺 区 的 变化 
幅度 并 不 明显 ,此 外 ,敏感 性 分 析 中 表 矿 比 辐射 率 
的 变化 对 消融 结果 的 影响 不 大 。 因 此 ,在 对 表 研 枚 
盖 型 冰川 进行 消融 模拟 时 ,应 高 度 关注 表 奔 导 热 系 
数 的 实地 观测 值 ,以 及 表 厂 表面 粗糙 度 和 表 太 反照 
率 的 取 值 ,这 些 表 奔 参数 的 微小 变动 会 对 消融 结 

产生 较 大 变化 。 由 于 本 文 未 对 所 有 表 磺 参数 开展 
实地 测量 ,其 率 定 值 可 能 与 青 冰 滩 72 号 冰川 的 实际 


数值 存在 一 定 的 差距 ,未 来 应 加 强 对 表 矿 属性 的 实 
地 观测 ,以 确定 其 时 空 变 化 和 不 确定 性 程度 。 
4.2. TARARE SY Wk JI kE 

为 进一步 理解 青 冰 滩 72 号 冰川 消融 区 域 的 能 
量 平衡 特征 ,本文 对 比分 析 了 不 同 区 域 山地 冰川 的 
能 量 和 消融 状况 ( 表 4)。 西 琼 台 兰 冰 川 和 科 其 喀 尔 
冰川 均 位 于 托 木 尔 峰 地 区 ,它们 与 青 冰 滩 72 号 冰川 
有 着 较为 相似 的 特点 ,冰川 面积 巨大 , 冰 崖 和 冰 面 
湖广 泛 发 育 。 西 琼 台 兰 冰川 面积 超过 100 km, vk 
面 差异 化 消融 现象 显著 ,1978 年 在 海拔 4000 m 处 表 
矿区 观测 的 能 量 收 入 来 源 中 , 净 辐 射 收 入 占 比 为 
73% , 感 热 通 量 占 比 23% ,凝结 潜 热 仅 占 4% 左 右 , 平 
均 消 融 能 量 为 77 Wem, M FEHI KR YKJ 
4200 mm 处 的 观测 站 显示 5 一 9 月 的 平均 气温 在 1%C 
左右 ,而 极端 最 高 气温 达到 了 12 % 左 右 , 太 阳 辐 射 强 
烈 , 有 利于 冰川 消融 ,与 此 同时 5 一 9 月 集中 了 全 年 
75% 的 降水 ,降低 了 冰川 能 量 交换 的 水 平 ",2004 年 
在 该 冰川 表 磺 区 模拟 的 消融 能 量 为 10.5 Wem, 
HER AS PA POSE BE Fa KC JT]? 的 能 量 通 量 分 布 状况 与 地 处 
天 山 的 青 冰 滩 72 号 冰川 较为 相似 , 仅 有 净 短 波 辐射 
和 洪 热 通 量 存 在 一 定 差异 ,可 能 原因 是 嘎 隆 拉 冰 川 
纬度 较 低 接收 了 更 多 的 短波 辐射 ,同时 受到 东亚 季 
风 环 流 的 影响 冰 面 水 汽 含量 较 高 从 而 产生 大 量 的 
洪 热 交换 。 与 青 冰 滩 72 号 冰川 相 比 ,阿尔 卑 斯 山脉 
的 Miage 冰川 接收 了 更 多 的 太阳 辐射 ,而 净 长 波 辐 
射 的 亏损 比 72 号 冰川 高 出 90% 以 上 ,传导 热 通 量 的 
能 量 消耗 也 十 分 巨大 ,最 终 导致 该 冰川 可 供 消 融 的 
能 量 不 足 5 Wm””。 喜 马 拉 雅 山脉 的 Lirung 冰 川 
纬度 低 于 30°N ,冰川 来 端 海拔 接近 4000 m,2013 年 
5 月 该 冰川 可 供 消融 的 能 量 为 123 W m? ,消融 程 
度 远 高 于 其 他 区 域 的 冰川 。 喀 喇 昆仑 山 所 处 纬度 
较 高 ,因此 该 地 冰川 接收 的 短波 辐射 显著 低 于 其 他 
区 域 , 除 此 之 外 ,Collier 等 ”在 喀 喇 昆仑 山 模拟 的 感 
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表 4 不 同 区 域 表 奢 履 盖 型 冰川 的 能 量 分 布 状 况 


Tab.4 Surface heat fluxes components of debris-covered glaciers in different regions 


青 冰 滩 72 号 冰川 FORGO BERR KSI] ” 呀 隆 拉 冰 川 Miage 冰 川 Lirung 冰 川 喀 喇 昆仑 山 
研究 尺度 单 点 尺度 区 域 尺度 
经 度 79*54'E 80*17'E 80*10'E 95°43'E 6°52'E 85°33'E 74°~77°E 
纬度 41°45'N 41°45'N 41°48'N 29°45'N 45°47'N 28°13'N 34°~37°N 
海拔 /m 3560-5986 3390-5800 3060-6342 3880-4800 1775-4800 4000-7132 3000-8000 
面积 km 5.74 107.1 83.6 2.7 11 6.3 18000 
表 矿 厚度 /m 0.12 一 1.8 0.25 0.23 一 0.01~1.07 
SIW +m?) 185.2 201 219 176 131.8 
LW +m?) -36.6 E vs -252 -70 -63.3 -85.6 
HI(W *m?) -92.1 -15 -43.5 -87.1 -75 11.3 4.3 
LEW *m?) -8.6 -54.7 -44.0 -25 -14 -1 -13.6 
GI(W +m”) - - - - -57 - 3 
PI(W *m?) -0.1 = = 0.1 x E 0 
MI(W * m?) 47.8 TI 10.5 63.8 3 123 39.9 
文献 来 源 本 文 [38-39] [41] [12] [28] [42] [43] 


E 
m o 


热 通 量 和 传导 热 通 量 均 为 正 值 ,而 其 他 单 点 太 度 的 
模拟 研究 则 显示 除 净 短波 辐射 外 ,其 余 能 量 一 般 为 
能 量 文 出 状态 ,这 表明 模型 尺度 和 空间 分 辨 率 的 差 
异 对 辐射 通 量 和 消融 模拟 结果 具有 重要 的 影响 。 
总 体 来 看 ,不 同 区 域 的 冰川 消融 状况 主要 取决 于 当 
地 的 环流 特征 与 地 表 气 象 条 件 ,同样 ,不同 的 能 量 
平衡 模型 以 及 气象 要 素 参 数 化 方案 的 差异 也 会 影 
响 模拟 结果 的 准确 性 。 


s ”结论 与 展望 


本 文 以 天 山 托 木 尔 峰 青 冰 滩 72 号 冰川 为 例 进 
行 消融 模拟 研究 ,该 冰川 所 在 区 域 草 藏 着 丰富 的 冰 
雪 水 资源 ,平均 每 年 产 流 约 63.4x10* m , Hor 24] 5096 
为 冰雪 融 水 , 旦 分布 着 众多 的 表 研 履 盖 型 冰川 ,这 
类 冰川 由 于 表 磺 覆盖 形成 差异 化 消融 ,对 冰川 物质 
平衡 和 径流 评估 带 来 一 定 的 挑战 ,因此 ,本 文 基于 
表 研 区 自动 气象 站 小 时 尺度 数据 作为 输入 ,采用 表 
磺 覆 盖 型 冰川 能 量 平衡 模型 开展 下 禾 冰 消融 模拟 
研究 ,结合 野外 实测 消融 和 表 磺 温度 数据 进行 验 
证 ,得 出 以 下 结论 : 

(1) 研究 期 内 能 量 平衡 模型 模拟 的 消融 量 为 
0.39 m w.e., 基 于 花 杆 实测 的 消融 量 为 0.40 m w.e., 
R=0.92, RMSE=+0.03 m w.e., 能 量 平衡 模型 取得 了 


注 :5 为 净 短波 辐射 ;L 为 净 长 波 辐射 ;为 感 热 通 量 ;LE 为 潜 热 通 量 ;G 为 传导 热 通 量 ;P 为 降水 带 来 的 热 通 量 ;MM 为 冰川 消融 能 量 ; 单 位 均 为 W: 


较 好 的 模拟 效果 ,模型 计算 的 表 研 表 面 温度 和 表 磺 
层 内 温度 也 与 野外 实测 温度 数据 具有 较 好 的 一 
致 性 。 

(2) 在 冰川 表 厂 区 的 能 量 收 支 过 程 中 , 奖 短 波 
辐射 是 唯一 的 能 量 收入 项 ,是 冰 面 的 主要 能 量 来 
源 , 其 他 能 量 项 均 处 于 亏损 状态 , 感 热 通 量 是 最 大 
的 能 量 文 出 项 (49.7% ) ,其 次 分 别 为 传导 热 通 量 ( 消 
融 耗 热 )(25.8% ) , fI UE f 5d (19.8% ) 和 潜 热 通 量 
(4.6% ) ,降水 热量 不 足 1% ,可 忽略 其 对 冰 面 能 量 交 
换 的 影响 。 

(3) 不 同 天 气 条 件 对 表 磺 区 气象 和 能 量变 化 有 
显著 影响 ,但 不 改变 它们 的 时 间 分 布 特征 ,入 射 短 
波 辐射 下 行 长 波 辐射 和 相对 湿度 受 云 量 的 影响 程 
度 较 大 , 阴 天 的 云 量 显著 削弱 净 短 波 辐射 通 量 ,其 
他 辐射 通 量 也 较 晴 天 减少 50% 左 右 ,了 明天 的 平均 消 
融 量 比 晴 天 减少 了 12%, 需 注意 云 量 对 冰川 能 量 收 
支 状 况 的 影响 。 

(4) 通过 对 表 厂 属性 进行 敏感 性 试验 ,模拟 的 
消融 量 对 导热 系数 的 变化 最 为 敏感 , 表 奢 反照 率 和 
表面 粗糙 度 的 变化 量 也 不 容 忽视 ,而 表 磺 比 辐射 率 
对 消融 模拟 影响 不 大 。 

本 文 仅 涉及 冰 面 单 点 的 消融 模拟 ,由 于 冰川 不 
同 区 域 的 气象 条 件 及 地 形 特 征 具 有 很 大 的 差异 性 ， 
为 了 厘清 整 条 冰川 乃至 流域 尺度 的 表 磺 覆盖 型 冰 
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川 的 消融 状况 ,未 来 应 从 以 下 三 个 方面 加 强大 区 域 
的 能 量 消融 和 动力 模拟 :(1) 完善 对 冰 面 气象 要 素 
的 时 空 变化 特征 研究 ,以 确定 气象 要 素 的 参数 化 方 
案 在 大 范围 区 域内 取得 最 佳 的 全 局 效果 ;(2) 加 强 
对 表 矿 厚度 以 及 表态 导热 系数 的 实地 观测 和 模拟 ， 
以 精确 驱动 消融 模型 ;(3) 构建 冰川 动力 学 参数 化 
案 ,并 耦合 能 量 平衡 模型 ,更 为 细致 地 研究 冰川 
动力 和 运动 特征 对 冰川 消融 的 影响 
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Simulation study of summer ablation in the debris area of Qingbingtan Glacier 
No. 72 in Mt. Tomor 
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Abstract: Debris-covered glaciers are widely distributed in Western China. Their ablation areas are covered by 
varying degrees of rock debris, and consequently, their melting statuses differ greatly when compared to debris- 
free glaciers. There is currently a need for melting simulations to better understand debris- covered glaciers. In 
this paper, driven by field meteorological data, an energy balance model for debris- covered glaciers has been 
used to simulate the energy and ablation in debris- covered areas of Qingbingtan Glacier No. 72 in Mt. Tomor, 
Tianshan. Based on the heat conduction process and the energy balance equation, the model calculates the debris 
surface temperature and the internal temperature of the debris, then estimates the subdebris melt using the internal 
debris temperature. The results showed that the modeled ablation was 0.39 m w.e. in the summer of 2008, and the 
simulation accuracy (R^ = 0.92, RMSE = +0.03 m w.e.) was higher when compared with the field data. The 
simulated debris temperatures at the surface and a depth of 10 cm inside the debris were also found to fit well 
with the measured data (R^ — 0.91 and 0.60, respectively). During energy exchange in the debris area, net 
shortwave radiation was the only energy income item, and sensible heat flux was the largest energy expenditure 
item (49.7%), followed by the heat conduction flux (ablation heat consumption) (25.8%), net longwave radiation 
(19.8%), and latent heat flux (4.6%), while precipitation heat was «196. Cloud cover had a significant impact on 
the meteorological and energy characteristics of the debris area. Under overcast conditions, the incoming 
shortwave radiation in the debris area decreased from 854 W - m? on sunny days to 587 W * m ^, while the 
downward longwave radiation and relative humidity increased, and the average ablation decreased by 12%, when 
compared with sunny days. In addition, the sensitivity analysis of the key parameters for debris shows that the 
simulated ablation is most sensitive to the changes in thermal conductivity, and the changes in albedo and surface 
roughness cannot be ignored. 
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